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Abstract 



We study neutrino physics in the context of a supersymmetric model where a continuous 
It-symmetry is identified with the total Lepton Number and one sneutrino can thus play the 



role of the down type Higgs. We show that R-brcaking effects communicated to the visible 
sector by Anomaly Mediation can reproduce neutrino masses and mixing solely via radiative 
contributions, without requiring any additional degree of freedom. In particular, a relatively 
large reactor angle (as recently observed by the Daya Bay collaboration) can be accommodated 
in ample regions of the parameter space. On the contrary, if the R-breaking is communicated to 
the visible sector by gravitational effects at the Planck scale, additional particles are necessary 
Qh to accommodate neutrino data. 
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1 Introduction 

Having already collected an integrated luminosity of 5 fb -1 , the LHC is starting to probe the 
nature of the (possible) UV completion of the Standard Model (SM). Supersymmetry (SUSY) is 
surely one of the best motivated SM extensions, since it elegantly solves the hierarchy problem. In 
the construction of the supersymmetric version of the SM (SSM), one finds dangerous operators 
that allow for proton decay. In order to forbid such operators, the common assumption is to 
enlarge the symmetry group to SU(2) x U(i) x G, where invariance under G forbids proton decay. 
Typically one assumes that G is a discrete group (R-parity R p ) under which ordinary particles are 
even while supersymmetric particles are odd. Beside forbidding Baryon (B) and Lepton (L) number 
violating operators that generate proton decay (and other flavor changing processes), an immediate 
consequence of R-parity is to make the Lightest Supersymmetric Particle (LSP) absolutely stable. 
Of course there are alternatives to R-parity (e.g. one can impose invariance under other discrete 
groups, under L and/or B, or one can extend R-parity to a continuous U(1)r [ll[2]), and one can 
even assume that proton decay is not forbidden, as in R-parity violating (RPV) theories (see [3] for 
a comprehensive review), where however the coefficients of the L and B violating operators must 
be strongly suppressed. 

The case of G = U(1)r, the continuous group that contains R-parity as Z2 subgroup, requires 
to go beyond the minimal scenario. Indeed, the R-symmetry forbids Majorana gaugino masses, 
but Dirac gaugino masses are allowed if the gauge sector of the theory is enlarged to the one of 
N = 2 SUSY, e.g. including adjoint superfields 3>f for each gauge group G{. If SUSY breaking 
is transmitted to the visible sector through a spurion D-term, (W„) = D'9 a , then a lagrangian 
term of the form -j-j J d 2 9(W'W^)^f generates Dirac mass terms of order nid ~ jr. R-symmetric 
models (#[6] represent an interesting possibility to explore for several reasons. First of all, gaugino 
one loop contributions to squared soft masses are finite [7j, so that the fine-tuning issue for the 
gluino is softened (see e.g. J8]). In addition, the LHC phenomenology is non standard, both due 
to the Dirac nature of the gluino [9-12 and to the presence of additional particles that can be 



rather easily detected [13] . Moreover, the Flavor Problem is also softened, since unsuppressed 
flavor changing terms are now allowed for sufficiently heavy gaugino Dirac masses |4||14|. 

As it has been recently explored in [61 and [17 1 rl it is not necessary to define the R-symmetry 
as the continuous symmetry containing R-parity. Indeed, proton stability could be ensured also 



identifying the R-symmetry with Lepton number [6] or with Baryon number 17 . Both scenarios 
violate R p , but proton stability is guaranteed without any suppressed coupling, since the model 
posses either an accidental standard Baryon or standard Lepton number. We will focus here on a 
scenario where the R-symmetry is identified with Lepton number. One of the distinctive feature of 
this idea is that it allows for a sneutrino to play the role of down- type Higgs |6rl Assuming the 
R-symmetry not to be spontaneously broken by the sneutrino vev, one is then forced to require 
vanishing Lepton number for the slepton doublet. This immediately implies that neutrino Majorana 
masses are forbidden by the R-symmetry and that the sneutrino vev, being unrelated to neutrino 
masses, can be large enough to give mass to the bottom quark. 

However, the R symmetry is not an exact symmetry, since an irreducible source of R- breaking 
(Xfrom now on) is given by the gravitino mass necessary to cancel the cosmological constant. This 
suggests a tight connection between neutrino physics and SUSY breaking. 



See also Il5j and 16 for earlier attempts. 



2 The idea of having a non zero sneutrino vev has been extensively explored in the literature, see 18 for examples. 
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Table 1: R-charge assignment for the chiral supermultiplets in our model. 

In the specific model presented in [6] , the R-symmetry was identified with the lepton number of 
a specific flavor, and only one non zero neutrino mass was generated. We want here to enlarge the 
R-symmetry to the total Lepton number to see whether this more realistic scenario can reproduce 
neutrino physics, analyzing in detail the parameter space compatible with the present experimental 
neutrino data. 



2 ^-symmetry as global lepton number 

Let us now describe our framework. We generalize the model of p] in such a way that the R- 
symmetry is identified with the global lepton number, XJ(1) R = U(1)l- In particular, all the 
R-charges of Lepton doublets and singlets are respectively fixed to and 2, see Table [2j The Rd 
electroweak doublet with R-charge 2, introduced to have an anomaly free framework, will play the 
role of an inert doublet (since we do not want the R-symmetry to be spontaneously broken), while 
the role of the usual down-type Higgs doublet will be played by a combination of sleptons, as we 
will explain later on. Since R-symmetry invariance is incompatible with Majorana gaugino masses, 
it is necessary to introduce three adjoint superfields, ^^ g -, that couple to the ordinary gauginos 
via D-term SUSY breaking u\ to generate Dirac masses. 

The most general superpotential compatible with the given R-charge assignment is: 



W 



I^H u Rd + H u QYjjU c + 2_^ KjkLiLjE^, + 



<jL 



+ Y^ KjkUQjD c k + \sH u §zR d + \ T H u <f> w R d . 

where Xijk = —Xjik from the antisymmetry of LiLj. 

The R-conserving SUSY breaking soft lagrangian is instead: 

£r = m 2 Hu h ] u h u + Y^ {m 2 Lij l\Ij + m 2 Rij e\e 5 ) - ^ b^hJi + 

ij i 

+ Yl ( m 9ij9iQj + m dijd\dj + mlijuluj) + M B Bi> § + M^tr(W^) . 



(2.1) 



(2.2) 



The R-symmetry cannot be an exact symmetry, since it is broken at least by the gravitino mass 
necessary to cancel the cosmological constant. To write down the ^ soft SUSY breaking lagrangian, 
we need an ansatz on how the R-breaking is communicated to the visible sector. A minimal 



scenario is to assume that gravity conserves the R-symmetry 19 , so that R-breaking effects are 
communicated to the visible sector only through Anomaly Mediation; however, we can also imagine 
that gravity effects at the Planck scale can break the R-symmetry. 

In the first case, which we will call Anomaly Mediation R-Breaking (AMRB) scenario, the soft 
R-breaking lagrangian is given by: 

C'jf = -CMajorana + C-A. + B^h u r d (2-3) 

where 



C A = AU^et + AfJ^dt + huqAVtf. (2.4) 



^Majorana ~ m B BB + m^ti(WW) + m g tr(gg) , 
1 A PPp c 4- A D 

The first term contains gaugino Majorana masses of order m — j^^ 3 l while the second one contains 
trilinear scalar interactions proportional to the supersymmetric Yukawa couplings. 

Turning to the case in which gravitational effects at the Planck scale break the R-symmetry 
(which we will call Planck Mediated R-Breaking (PMRB) scenario), the R-breaking structure is 
much richer than in the previous case, since now all the operators suppressed by some power of the 
Planck scale can contribute. The R-conserving superpotential and soft SUSY breaking lagrangian, 



Eqs.(2.1), (2.2), are corrected by the following ^-contributions: 



i 
£ W- = ^Majorana + ^A + BfJluTd- (2.5) 



The Ji soft SUSY breaking contribution has the same structure as in Eq. (2.4), but now we simply 
expect all the terms generated to be of order of the gravitino mass 77J3/2 and the A-terms not to 
be aligned to the supersymmetric Yukawa couplings. Let us notice the appearance of //-terms and 
Majorana masses for the Adjoint Fermions, also of order 7713/2- As we will see, they will play an 
essential role in neutrino physics. 

Let us now study how electroweak symmetry breaking works in this framework and how fermions 



get masses. Since all the sleptons have a b^ term, Eq. (2.2), in a general basis all sneutrinos will 
get a vev. However, we can use the freedom to rotate slepton fields to work in a "single vev basis" 
where just one sneutrino gets a vevQ We will denote with A, B, C the flavor indexes in this basis, 
with La referring to the doublet that plays the role of the down-type Higgs. The superpotential 
can be rewritten as: 

W = ^H u R d + \ s H u <f> B R d + \ T H u $ w R d + WVukawa + TU tri i inear , (2.6) 

with 

Yukawa = VbLaL b E% + y c L A LcE c c + yfL A Q t Df + H U QY U U C , 

Wti-ilinear = 2-*/ ^BCjLbLcEj + ^ j (XsjjLbQjD^ + XgjLcQiDj), (2.7) 

i=A,B,C ij 



20 



3 See 
4 This is 



for the exact expressions of gaugino Majorana masses and A-terms in Anomaly Mediation, 
similar to what happens in RPV SUSY models, 31. 



where y B = A abb and yc = Aacc- in the new basis, the R-conserving soft lagrangian of Eq. (2.2) 



maintains the same form, while the R-breaking ones of Eqs. (2.4), (2.5) now read 

C-A = £lR + Arilinean (2-8) 

with 

£lr = A B l A £ B e c B + A c i A I c ~e c c + A bc I A Ib^c + A CB iAice c B Af ) £ A q i d c i + h u qA u uf, 

Arilincar = A BCi £ B £ C e\ + A^£ B q4) + A' Ci f£ C qUd). (2.9) 

The first term gives slepton and squark left-right mixing qJ while the second term contains trilinear 
scalar interactions that do not involve the slepton that takes vev. Let us stress that the gaugino 
Majorana masses and the scalar left/right mixing will play a crucial role in the generation of 
neutrino masses. 

The analysis of the scalar potential can be done along the line of Ref. [3] , although in our case 
the situation is more involved. Indeed, when the left handed slepton soft squared mass matrix is 
not flavor universal, a mixing between the sneutrino that takes vev and the other two is in principle 
possible, so that we expect the physical Higgs to be an admixture of all the three sneutrinos. On 
the contrary, when the squared mass matrix is flavor universal, the resulting scalar potential is the 
usual one [6|. We assume here for simplicity that, at leading order, the soft squared mass matrix 
is flavor universal, deferring to a future work the analysis of the non flavor universal case. 



From Eq. (2.7) it is immediate to notice that the charged lepton of flavor A cannot acquire mass 
trough a SUSY invariant Yukawa term as the operator £ A £ A e A ^ s nun due *° * ne SU(2) invariance. 
Therefore, a mass for the lepton £ A must be generated by a hard SUSY breaking sector through 
couplings between messengers and leptonic superfields J6j. However, in the present scenario, this 
sector will generate hard Yukawa couplings also for the B and C flavors. 

If we assume that the main contribution to £ B) c masses comes from the supersymmetric Yukawa 
couplings, the additional contribution from the hard sector must somehow be suppressed. This 
makes A = e the simplest possibility. Indeed, if A = r, the r lepton mass must be generated by 
the hard sector, while the hard contribution to the other masses must be suppressed (for example 
requiring the hard Yukawa couplings yij to satisfy yij <C 10~ 6 ). This corresponds to assuming a 
large hierarchy between the hard Yukawa couplings. The same line of reasoning can be applied in 
the A = fi case. If instead A = e, a hard Yukawa contribution which generates Yukawa couplings 
of order y e ~ 0(1O -6 ) for all the charged leptons does not give a too large contribution to the fi 
and r masses, while providing the correct order of magnitude for an electron mass. Since in this 
case there is no need to introduce any large hierarchy in the new sector, it appears a more natural 
choice. A possible example of hard Yukawa sector is given in |6|; however, let us stress that, since 
we will left largely undetermined this sector, in what follows we will analyze also the cases in which 
A^e. 

As a last comment, let us stress that the interaction terms of Patrilinear (which are not present 
in (6l) closely resemble the trilinear interaction terms that appear in RPV theories (31. However, in 
our case all the off diagonal terms involving the flavor A, AV ■, are zero in the single vev basis, so 

that the number of parameters is reduced. Moreover, coupling of the type X Aii now play the role 
of Yukawa couplings and are not free parameters. We conclude that our scenario is a variation of 



3 In the AMRB scenario the off diagonal terms Abc, Acb are zero. 



RPV models (with less parameters) , although as we will see a larger amount of R-parity violation 
in the neutrino sector than in the standard case will be allowed. 

2.1 Electroweak precision measurements and flavor constraints 

Let us now discuss in turn the experimental constraints coming from Electroweak Precision Mea- 
surements (EWPM) and from flavor physics. 

One of the distinctive feature of models where the R-symmetry is identified with Lepton Number, 
is that all the supersymmetric partners, with the exception of charged sleptons and sneutrinos, have 
a non vanishing lepton number. As a consequence, charged leptons and neutrinos can mix with 
the "new" spin 1/2 leptons (Dirac gauginos and higgsinos). A priori, the neutralino mass matrix 
is a 9 x 9 squared matrix, while the chargino mass matrix is a 12 x 12 square matrix. However, in 
the single vev basis, the leptons of flavors B and C do not mix with any other fermion, so that the 
effective matrix is the same as in I6J, to which we refer for a detailed analysis of the mass eigenstates. 
The important point to stress for our purpose is that in the R-symmetric limit all neutrinos are 
massless. Also, the same bounds on the sneutrino vev coming from the bounds on the coupling of 
the Z boson to charged leptons apply, i.e. for M^ ~ 1 TeV one should have va i$ 40 GeV. 



Let us now turn to the bounds on trilinear couplings appearing in WVukawa an d Patrilinear 21 
Since these are RPV couplings, we refer to [3j[2l] for a detailed description of the origin of the 
various bounds. It is interesting to notice that our framework has distinctive differences both with 
the model of p] and with the standard RPV SUSY. 

On the one hand, Lepton Flavor Violating (LFV) processes are allowed in our framework but 
not in p] , and the same is true also for semileptonic meson decays (such as rare decays of B and K 
mesons), unless we assume alignment between the matrices (A^^)^- and the quark mass matrix. 

On the other hand, even though our situation is more similar to the standard RPV SUSY, 
some bounds have a different interpretation. In particular, bounds that involve a product between 
two trilinear couplings can now involve one Yukawa coupling. In order to maximise the parameter 
space for the sneutrino vev we read these bounds as vev dependent constraints on the trilinear 
couplings appearing in Patrilinear- For example, when A = e, the LFV process \x — > cy puts a 

a ( rn i \ ^ 
bound | A233A133I < 2.3 x 10 ( 100 jjgy J (assuming for simplicity degenerate slepton masses). 

Since in our model A133 = ^ is the r Yukawa coupling, the bound can be restated as | A233I < 

100 ^ eV j ~ 0.002 - 0.07 for m^ ~ 200 GeV and v e = (10 - 80) GeV. At the 
same time, it is true that, among the constraints that involve only one trilinear coupling and not a 
product, some will refer to bounds on Yukawa couplings, implying thus a bound on the sneutrino 
vev. 

In the following we will always assume that trilinear couplings not directly related to neutrino 
physics are always small enough to satisfy all the experimental constraints. 

3 Neutrino physics and U(1)r lepton number 

In our model the R-symmetry is identified with the global Lepton number, so that U(1)r breaking 
corresponds to Lepton Number breaking. In the following section we will discuss how neutrino 
masses and mixing are generated from R-symmetry breaking effects. The problem of neutrino 
masses in models with an R-symmetry have been studied both for Majorana [22] and Dirac [23] 
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Table 2: Fits to neutrino oscillation data. Where two different values are present 
for one parameter, upper and lower row refer respectively to Normal and Inverted 
Hierarchy. 



neutrinos. Both scenarios require to enlarge the particle content of the model introducing right 
handed neutrinos. Indeed, in the standard R-symmetric scenario (4^5], there is no natural connec- 
tion between the R-breaking and Majorana neutrino masses, since these are allowed by R-symmetry 
(all lepton superfields have R-charge 1). A priori, however, R-symmetry does not forbid Dirac 
masses either, since their presence depend on the R-charge assignment of right-handed neutrinos. 
This makes the connection between neutrino Dirac masses and R-symmetry breaking less stringent. 

On the contrary, in our scenario there is a clear connection between Majorana neutrino masses 
and R-breaking effects, since such Majorana masses are clearly incompatible with the U(1)r sym- 
metry. In this way, in principle we don't need to introduce any additional particle (i.e. right-handed 
neutrinos) in order to generate non zero masses. While, as we will see, this will be true for AMRB, 
in the case of PMRB additional structure will be necessary in order to reproduce neutrino masses 
and mixing, making this scenario less compelling. 

Let us stress again that in our scenario the scale at which Lepton Number is broken is deeply 
connected with the scale of supersymmetry breaking through the gravitino mass, while in general 
the Majorana neutrino masses generated through the Weinberg operator call for a very large scale, 
which may or may not be connected to the scale of supersymmetry breaking. 



3.1 Neutrino masses and mixings 

Before analyzing the neutrino phenomenology in our framework, let us briefly summarize some 
features of a general neutrino mass matrix. 

As it is well known, the neutrino mass matrix is largely undetermined, since we lack of informa- 
tion on the absolute neutrino mass scale and on the hierarchy between the mass eigenstates. For 
three active neutrinos, the present data are summarized in Table [2] 

At the same time, CMB data point towards Yli m i>,i ~ 0.6 eV (see e.g. |27j), from which one 
can infer a loose upper bound mughtest ^S 0.1 eV for both hierarchies. Using data in the expression 
of the neutrino mass matrix in terms of masses and mixing, we expect the following general form 
for the mass matrix (in the (v e , v^^v^) basis): 



(*) Normal Hierarchy: 



M v ~ m s J naU I 0(e) 0(1) 0(1) 1 , Af„ ~ m^ e I 0(e) 0(1) O(^) ) ; (3.1) 





0(1) / 



(*) Inverted Hierarchy: 



/0(1) 0(e) 0(e)\ /0(1) 0(e) 0(e) \ 

M u ~ mr" 0(e) 0(1) 0(1) , Af„ ~ m^ e 0(e) 0(1) 0{y/t) ; (3.2) 
\0(e) 0(1) 0(1)J \0(e) O(^e) 0(1) J 

The matrices on the left and on the right refer respectively to a small (of order 0(1O -5 eV)) and 
large (of order 0(1O _1 eV)) lightest neutrino mass. The exact value of the coefficients depends on 
the chosen values of neutrino masses and mixing angles (we discard here the dependence on Dirac 
and Majorana phases); however, typical order of magnitudes for the elements are 

m small _ O ( 10 -2 eV ) ; m^ 6 ~ ©(HT 1 eV), e~O(10~ 1 eV), e ~ 0{l^ eV) 

In what follows, taking the approach of [28], we will focus on specific forms for the neutrino mass 
matrix that we consider representative of the different phenomenological scenarios. In particular, 
we will focus on the two following matrices, representative respectively of the Normal and Inverted 
Hierarchy cases for small lightest neutrino mass: 
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eV ; mjf ~ 10" z -0.54 2.19 -2.36 eV . (3.3) 
\-0.52 -2.36 2.8 / 

We do not show here the corresponding matrices for the large lightest neutrino mass scenario 
because, as we will explain later on, they can be reproduced only in a very small region of parameter 
space. To construct the previous matrices, we fixed the lightest neutrino mass to 2 x 10~ 5 eV, 
while the other parameters are fixed as follows: A?nf 2 ~ 7.6 x 10~ 5 eV 2 , Am 2 3 ~ 2.4 x 10~ 3 eV 2 , 
sin 2 #i2 ^ 0.3, sin 2 #23 — 0.47, sin 2 #13 ~ 0.024, i.e. we take #13 ~ 9° as recently observed by the 
Daya Bay collaboration |29 . For simplicity, we have also assumed a vanishing CP violating phase. 



3.2 Neutrino physics in AMRB 



Inspecting Eq. ( |2.3| ), it is clear that the gaugino Majorana masses contribute to the neutralino- 
neutrino mass matrix. This resembles what happens in RPV theories with bilinear terms [3] , where 
one neutrino gets a non zero mass already at tree level through its mixing with gauginos. 

On the contrary, in this scenario all neutrinos remain massless at tree level. This is a striking 
difference with respect to the RPV case, and can be understood considering the approximate 
eigenstates of the neutrino mass matrix (calculated e.g. using the usual seesaw formula): 
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Figure 1: Allowed region (colored) in parameter space for the flavor assignment A = e, 
B = |U and C = r in the case of Inverted Hierarchy. 



The B and C flavors are by themselves approximate eigenstates and cannot get mass through 
a mixing with gauginos. At the same time, the flavor A mixes only with Higgsinos and adjoint 
fermions, so that the absence of mixing with gauginos and of Majorana masses for the adjoint 
fermions prevents va from getting a tree level mass. 

It is now clear that, in the AMRB scenario, the only possibility for neutrinos to acquire a mass 
is through loop effects. In the (1^4, vb, vc) basis, the main contributions at 1-loop are given by J3J: 

• Loops with two supersymmetric trilinear couplings and one mass insertion in the scalar prop- 
agator due to Anomaly Mediation. 

Since this term is proportional to the mass of the fermion circulating in the loop, the dominant 
contributions are given by bottom quark, strange quark and tau lepton [j 



6 We neglect the muon contribution because, due to color factors, it is subdominant with respect to the strange 
quark contribution. 



They are given by: 
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where \ Aii = (md)i/vA is the i down-quark Yukawa coupling, m~ h is the common left handed 



and right handed sbottom mass scale, fib is the bottom /3-function 20 , and for simplicity we 
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In the lepton sector, the main contribution is given by 
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where A^33 = m T /vA is the tau Yukawa coupling and /3 T the tau /3-function 20 . 

Loops with two gauge couplings and one Majorana mass insertion in the gaugino propagator: 
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where M-^ is the Dirac Wino mass and we have used the Anomaly Mediation contribution 



to the Majorana Wino mass: m 
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In the previous equations we neglected the mixing of v a with the adjoint gauginos (see Eq. (3.4)): 
this is consistent in the portion of parameter space we will consider in the following numerical 
analysis. 

Barring special relationship between the parameters involved, the neutrino mass matrix has 
now three non zero eigenvalues. These depend on free parameters (trilinear RPV couplings and 
gravitino mass), that can be chosen to fit the experimental data, but also on gauge couplings and 
masses which are constrained by collider experimental bounds Q 

As already stressed, our scenario is a particular case of RPV SUSY (in particular the loop 
contributions are the same in both cases), so that it is interesting to compare the two situations. 
Usually in RPV scenarios the left /right sparticle mixing and the Majorana gauginos mass are at the 
EW scale, while in our case they are proportional to the gravitino mass and can be subleading for 
small supersymmetry breaking scale. This implies that while usually one needs to suppress too large 



7 In what follows we will always take as reference a "natural" spectrum for the supersymmetric partners, with 
only the squarks of the third generation below the TeV scale, while all other superparticle masses can be above the 
TeV scale. At the moment, the experimental bounds on this kind of spectrum is less severe than those obtained for 
almost degenerate squarks pi. Note that Dirac gauginos have an improved naturalness with respect to Majorana 
gauginos [?], and this allows us to have a natural gluino above the TeV scale and a heavier Wino. In what follows, 
we will take the Dirac Wino mass up to 10 TeV. 
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Table 3: Range of parameters used in the scan of Sees. 3.2.1||3.2^2 



loop contributions to neutrino masses putting severe upper bounds on the trilinear couplings J3J, 
in our case the upper bound is translated on the gravitino mass (with trilinear couplings usually 
allowed to saturate the bounds from EWPM and flavor physics, see Sec. |2.1[ ). 

A loose upper bound on the gravitino mass can be derived from cosmological considerations. 
Indeed, as already stressed, the absolute neutrino mass scale is bounded from above from CMB 
measurements, m v < 0.6 eV. This readily translates into an upper bound on the gravitino mass, 
which can be roughly estimate as follows. Since rriAA is the only entry in the neutrino mass mass 
matrix that do not depend on trilinear couplings, we can use it to roughly set the largest neutrino 



eigenvalue scale. For typical value of sparticle masses 



m~ h , < 1 TeV, 



M w < 10 TeV) we obtain 



in 



3/2 



< 0.5 GeV. 



We will now study in detail whether, in the AMRB scenario, the phenomenological neutrino mass 
matrices can be reproduced in the case where the flavor A is the either electron, muon or tau. 



3.2.1 A = e: Electronic Higgs 

In this case we assign A = e, B = u, C = r. We perform our numerical scan for the parameters 
of Table [3j requiring the other variables to reproduce the phenomenological matrices and imposing 
the constraints of [21] . For simplicity, we have assumed degeneracy between LH and RH sparticles, 
and a full family degeneracy in the slepton sector. Strictly speaking, this simplification implies 
that, barring accidental cancellations, a natural common slepton mass cannot be too large, since 
it enters in the determination of the Z mass through the minimization of the scalar potential. 
However, keeping in mind that only the LH slepton mass matrix affects the Higgs sector, and to 
have an idea of the general behavior of the model, we allow the common slepton mass to assume 
also larger values. 

The main result of this section is that while the Normal Hierarchy case can be reproduced 
only in a very small region of the parameter space (corresponding to v e ~ 100 GeV and rather 
large Dirac Wino masses, Myy > 5 TeV), a much larger portion of parameter space is available for 
Inverted Hierarchy. This can be understood looking at the phenomenological matrices of Eq. (3.3): 



the m ee entry in the Normal Hierarchy case is about one order of magnitude smaller than the one 
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of the Inverted Hierarchy case. For this to happen, one needs large sneutrino vev and large Wino 
mass. This can be seen noting that we can parametrize m ee as 



m ee oc I h /3v e J m 3/2 (3.8) 

where the first term comes from the squark and slepton loops (a ~ 1/m 2 ) while the second one is 
due to the Wino loop (/3 ~ 1/M^). A large vev can suppress the first term, while a large Wino 
mass can suppress the largeness of the second one. The available parameter space for Inverted 
Hierarchy is shown in Fig. [TJ where the allowed region is the colored one. We do not show plots 
on the sparticle and Wino mass planes since these parameters are practically unconstrained. As 
can be seen, in the squark sector the diagonal trilinear couplings A 333 , X' 233 are rather small, both 
at most of order 0(1CP 2 ), while the off-diagonal trilinear couplings A 332 , A 23 2 can be large, up to 
O(10~ 1 ). In the lepton sector we have again couplings A233, A231 at most of order O(10 _1 ). 

Another interesting consequence of our analysis is that we can set a more precise range on the 
gravitino mass, 

1 MeV < m 3/2 < 100 MeV, (3.9) 

Furthermore, we also have an indication on the sneutrino vev: we can fit the neutrino mass matrix 
in our framework only if the sneutrino vev is somewhat large, v e > 30 GeV, i.e. tan/3 = ^ < 6. 
Let us also notice that for larger sneutrino vev, a larger gravitino mass is allowed. This can be 



understood from Eq. (3.8), from which it is clear that for small sneutrino vevs the term between 
brackets can be large, so that in general a small gravitino mass is needed to suppress this entry. 
On the contrary, for larger values of the vev the term between brackets is more suppressed, and a 
larger gravitino mass is allowed. 

A comment on the situation for larger lightest mnght es t is in order. We have explicitly checked 
the situation for mn g htest — 0.1 eV, finding that only in very a small region of parameter space the 
phenomenological neutrino mass matrix can be reproduced. However, let us stress that in this case 
approximately the same region of parameter space can reproduce both Hierarchies, since now the 



typical form of the mass matrix in the two cases is similar (Eqs. (3.1), (3. 2)) 



3.2.2 A = /j,,t: Muon and Tau Higgs 

As pointed out in Sec. [2l we consider the case of an Electronic-Higgs (^4 = e) more motivated from 
the point of view of the generation of the hard Yukawa couplings. However, for completeness we 
study also other possibilities. In particular, as we will see, the A = [i case offers an interesting 
different phenomenological situation with respect to the A = e case. 

Let us start with A = fi, B = e and C = r. In this case Eq. (3.8) is valid for m w , which 



is similar for the two hierarchies. Thus, in general we expect that, unlike what happens in the 
A = e case, both Hierarchies should be reproduced. This is indeed what happens, as confirmed by 
the scan performed for the parameters of Table [3] within the same approximations described in the 
previous section. 
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Figure 2: Allowed region (colored) in parameter space for the flavor assignment A = fi, 
B = e and C = r in the case of Normal Hierarchy. 



The results are shown in Figs. [2]|3j Also in this case we have checked that increasing the lightest 
neutrino mass diminishes drastically the available parameter space (although also in this case both 
Hierarchies can be accommodated). 

As can be seen from the plots, the range of parameters is roughly the same as the A = e case, 
although some details can change. An exception is given by the lepton trilinear coupling A133, 
which is now allowed to be also of 0(1). Regarding the muon-sneutrino vev and the gravitino 
mass, interestingly the situation does not change much with respect to the A = e case: we conclude 



that the bounds of Eq. (3.9) are rather typical, for small neutrino masses, while they are no longer 



valid increasing the lightest neutrino mass. 

Let us now comment on the Tau-Higgs case, i.e. A = r. We have performed our analysis 
both in the approximation of vanishing (i) and non vanishing (ii) muon mass. In the case (i) 
there is no contribution from loops involving sleptons, so that one can solve for the Dirac Wino 
mass instead of scanning on it. The results of our scan show that a solution compatible with the 
phenomenological mass matrices, Eq. ( |3.3| ) requires either very large trilinear couplings or very 
large Wino masses (well above 100 TeV). While the first possibility is excluded by the bounds 
coming from EWPM [21], the second one is in principle viable. However, since as we already 
pointed out we want to stick to a spectrum which is not too unnatural, we consider this possibility 
at best marginal. 
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Figure 3: Allowed region (colored) in parameter space for the flavor assignment A = /j,, 
B = e and C = r in the case of Inverted Hierarchy. 

In the case (ii) there is a non vanishing slepton loop contribution, in such a way that the scan 
on parameter space is quite similar to those of the two previous Sections (with the exception that 
in this case one of the two trilinear coupling constants involved is the muon- Yukawa coupling, so 



that there is no need to scan over it). Nevertheless, also in this case compatibility with Eq. (3.3) 
requires trilinear couplings incompatible with the bounds of |21|. 



The situation is summarized in Table [3j and the conclusion is that the case A = t cannot reproduce 
neither a Normal nor a Inverted Hierarchy spectrum. 

3.3 Neutrino physics in PMRB 

Let us now turn to the case where gravitational effects also break the U(1)r symmetry rl The main 
difference with the previous case is that now two non zero neutrino masses are generated at tree 
level. To understand this, let us consider the mixing among fermions in the neutralino sector. In 



We thank T. Gregoire, P. Kumar and E. Ponton for discussions on this whole section. 
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the R-symmetric limit, the R = — 1 mass eigenstates are well approximated by: 



gvA , _ g'vA _,. , gvAVu^r 

IM W W V2M 

gvA . gv u 



VA * ^ " V2M^ + V^e^ + V^ 



^ * ^ ~ vaS^ 1 " 7^' (3 0) 

while i^b and i/c do not mix, as we have already noticed. The R = 1 states are instead: 

W'~W- ^ T h d , & ~B- / As h d . (3.11) 

The inclusion of R-breaking effects generates new mixing terms for all neutrinos: 

/ gvA , g'vA , . gvAVu^T Wr 

VA * VA - vm;** + v^b^ + 7^ - 7" " d ' 

^ ~ v B -—h d , (3.12) 

A 1 

which in turn produce a mass term for 1/4 and mixing terms mAB, m AC'- 

foV A v u ( g\ T g'Xs\ 

rriAA = V2 — — fi A , 

ix \M W M B J 

vav u ( ' gXr g'Xs\ 

mAB = 7vU¥~^J w " 

vav u ( gXr g'\s\ ,„ rg ^ 

mAC = ttA^-^)^ (3 ' 13) 

Furthermore, a Majorana mass for the adjoint gauginos is generated, and through it the neutrino 
va acquire an additional mass term: 

(QVAtf . ( g'vA, 2 (o -\A\ 

mA = i M^ )mT + { MB- )mS - (3 - 14) 

This is an example of inverse seesaw mechanism [30], where the role of the right handed Dirac 
neutrinos is played by the Dirac gauginos. Therefore, the tree level mass matrix in the PMRB 
scenario is: 

faV A v u (g\T g'Xs\ . (g'vA\ 2 . ( gvA\ 2 
mAA = ^^T (% - ft ) " A + life j ms+ UrJ "»■ 

vav u ( gX T g'\s\ 

mAB = Vfefe"*J" B ' 

vav u ( gX T g'Xs\ /„ rK x 

mAC = wA^w~^b)" c ' (315) 
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which has indeed just a zero eigenvalue. 

Let us first of all discuss the upper bound on the gravitino mass imposed by the condition 
m v ;$ 0.6 eV. Looking at the non zero entries of the mass matrix we see that in general the upper 
bound depends on the value of \t,s- As in the AMRB case, we focus on the tuaa entry. When 
the first term is negligible, the inverse seesaw term gives an upper bound 7713/2 < 1 — 10 keV for 
Mw ~ 1 TeV and va ;$ 100 GeV. On the other hand, when the first term cannot be neglected, 
it dominates over the term coming from the inverse seesaw, and the upper bound now reads 
< \ which can be more stringent than in the previous case (depending on the value 



m 



3/2 ^ -j, 



S,T 



of Xt,s)- We conclude that, under these assumptions, in PMRB the upper bound on the gravitino 
mass can be significantly lower than the one of the AMRB scenario. 

Let us now explain why in this case fitting neutrino physics calls for the introduction of a new 



sector in the model. Inspecting the phenomenological mass matrices of Eq. 3.3, we see that both 



hierarchies require leading order entries in the \x — t sector, which cannot be accommodated by the 



mass matrix (3.15). This is true for any choice of the flavor A. At the same time, we expect loop 
factors to be much smaller than the tree level entries, so that the overall picture cannot be modified 
too much. This calls for the introduction of a new sector in the model. We can wonder what is the 
minimal sector able generate neutrino masses and mixing. First of all we would like to generate 
neutrino physics without the need for a new source of R-breaking. This means we should consider 
a mechanism that generates neutrino masses and mixing when the lepton number is broken at very 
low scale (the keV gravitino mass). The minimal possibility we can think of is an inverse seesaw 
mechanism with additional electroweak singlets Q 

Therefore, we introduce a right handed Dirac neutrino (two singlets S and S with R = and R = 2 
respectively) and the following terms in the superpotential: 

W = Y^ hSH u Li + M S SS . (3.16) 

i 

Each singlet gets a Majorana mass of order of the gravitino mass trough R-breaking effects, 
and this generates a Majorana neutrinos mass of order m u ~ -^p-TOg. An interesting possibility 
for the Dirac mass Ms is the TeV scale, since this opens up a link between neutrino physics and 
LHC physics; however, a complete analysis of this situation is beyond the scope of the paper and 
we defer it to a future work. 

4 Conclusions 

With a luminosity of about 5 fb _1 already collected by the LHC, and without any hint of signal 
so far, the available parameter space of standard supersymmetric models is getting more and 
more constrained. This motivates the study of a larger portion of the weak scale supersymmetry 
landscape. Since neutrino physics can be a natural probe into new physics, it is natural to ask 
whether or not, given a specific framework, neutrino masses and mixing can be accommodated. In 
this work we have studied a supersymmetric scenario where a continuous R-symmetry is identified 
with the total Lepton Number, so that a possible connection to neutrino physics is immediate. 
In particular, we have found that neutrino physics is strongly connected with the mechanism of 
R-symmetry breaking, which in turn is related to supersymmetry breaking. 



3 In general, such Singlets may be present in the sector that generates the hard Yukawa coupling [6|. 
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When R-symmetry breaking effects are communicated to the visible sector solely via Anomaly 
Mediation, all neutrinos acquire mass at 1-loop level. The hierarchy that can be reproduced depends 
crucially on the flavor of the sneutrino that gets a vev and plays the role of down type Higgs. For 
small values of the lightest neutrino mass, and for A = e, the case of Normal Hierarchy is disfavored, 
since it can be reproduced only in a very limited portion of the parameter space. On the contrary, 
for A = fj,, both hierarchies can be fitted in a consistent portion of parameter space. Finally, 
for A = t, we are not able to reproduce neutrino phenomenology solely via loop effects. The 
situation changes increasing the lightest neutrino mass, since in this case both hierarchies can be 
accommodated for A = e and A = ^ (but not for A = r), but only in a limited region of parameter 
space. 

Another possibility is that R-breaking effects are communicated to the visible sector at the 
Planck scale. In this case two non vanishing neutrino masses are generated at tree level, but with 
a pattern that does not allow to reproduce the phenomenological matrices studied. Since loop 
effects give subdominant contributions and cannot change the overall picture, we conclude that 
a new sector must be added to the theory in order to reproduce neutrino physics. The minimal 
possibility is to introduce additional singlets (that can however be already present in the sector 
that generates the hard Yukawa couplings) in order to have an inverse seesaw mechanism. The 
study of this possibility is however beyond the scope of this paper. 

Since neutrino physics selects a particular region of the parameter space of the model, some 
consequences on Dark Matter and collider physics can be inferred. The cosmological upper bound on 
the sum of neutrino masses translates into an upper bound on the gravitino mass, m 3 / 2 < 0.5 GeV 
for AMRB (with a more precise range selected by the neutrino mass matrix fit, 1713/2 — 1 MeV — 
100 MeV)), and 1713/2 < 10 keV for PMRB. In both scenarios the gravitino lifetime is long enough 



to evade all experimental bounds, so that it can be a Dark Matter candidate 31 . 

Furthermore, neutrino physics selects also a preferred order of magnitude for the trilinear cou- 
plings both in the lepton and quark sector (with the general indication that the off diagonal cou- 
plings are larger that the diagonal ones). This can have important consequences for LHC physics. 
Indeed, one can expect squarks generation changing decays (as bj, — > vbsr or II — > ^sr) to dom- 
inate over the corresponding generation conserving decays (Sl —> vb^r or ti, — > e^bjt). A similar 
conclusion applies in the slepton sector, with decays like vr — > bs or £3 — > si generally dominating 



over i>b —> bb or e~B — > bt. We defer to a future work [31] the detailed analysis of possible signals. 
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